Predictors of psychological stress and behavioural diversity among captive red panda in Indian zoos and their implications for global captive management

At high elevations, biodiversity is at elevated risk from extinctions due to rapid environmental changes. In the most of its range in Himalayas, the red panda, an endangered species, is struggling to survive in the wild, and a global captive breeding programme has been launched to conserve the species. Because captivity can have negative impacts on animals, reducing the chance of successful reintroduction, we investigated the predictors of stereotyped behaviour and behavioural diversity of red panda (n = 26), and the effect of stereotypy on their behavioural diversity in three Indian zoos. Multivariate analysis showed that stereotypy increased with density of logs on the ground, age and higher among pandas in zoo 3 compared to zoo 2, but decreased with number of nests, sociality, tree density and tree height used. Similarly, behavioural diversity increased with log density, but decreased among pandas in zoo 2 compared to zoo 1, during summer compared to winter, and also with ambient temperature, stereotypy, tree density, and tree height used. The relationship between stereotypy and behavioural diversity was negative, but not significant. Provision of a greater density of trees was associated with higher levels of daytime inactivity. Findings from this study have global relevance, as the captive red panda experiences similar welfare issues around the world, and our data provide empirical support for some existing guidelines for red panda husbandry.


Results
Among 26 red pandas observed during daylight hours, 24 individuals showed stereotypic behaviours. Pacing, tongue flicking and position circling were the three behaviours observed. Position circling was observed among 14 individuals at Zoo 1 (54% of total assessed population across three zoos); this is the first report of this behaviour in pandas. The mean proportion of time spent in stereotyped behaviour was 2.0% ± 0.625% (54.5 min), with a mean behavioural diversity index of 1.07 ± 0.032 (using Shannon's H measure). Both stereotypy and behavioural diversity varied with some of the biological and environmental factors, and also differed between zoos. In the current study the overall activity budget of pandas was: inactive 59%, non-stereotyped activity 39% and stereotypy 2%.

Extent of stereotypic behaviour in relation to biological and environmental factors: univariate analysis.
Among the 15 independent factors tested against the extent of stereotypy, age, zoo and frequency of feed showed significant variation on univariate tests (Table 1). For example, pandas showed significant variation in stereotypy across the three different zoos (p < 0.05), with the lowest level in Zoo 2 followed by Zoo 1 and the highest in Zoo 3. Further, Scheffé's post-hoc test revealed that there was no significant difference in stereotypy between Zoo 1 and Zoo 2 (p > 0.05), while the individuals in Zoo 3 stereotyped more than that of Zoo 1 (p < 0.05) and Zoo 2 (p < 0.05). The extent of stereotypy was higher in pandas fed once compared to those fed twice (p < 0.05). Sub-adult pandas showed more stereotype than adults and cubs (p < 0.05). Stereotypy increased significantly with the density of logs on the ground (β = 27, R 2 = 0.32, F = 12.63, p < 0.05, Fig. 1); however, we found that the relation between log density and stereotypy might be deceptive, because four individuals that showed highest stereotypy were housed in small enclosures with no trees but more logs on the ground. Removing  Behavioural diversity in relation to biological and environmental factors: univariate analysis. Among the independent factors tested, behavioural diversity varied significantly only with zoo ( Table 1). The pandas managed in different zoos exhibited different levels of behavioural diversity (p < 0.05). Analysis using Scheffé's test showed that behavioural diversity differed significantly only between the individuals in Zoo 1 and Zoo 2 (p < 0.05). Univariate linear regression revealed that behavioural diversity decreased with Influence of biological and environmental factors on stereotypy and behavioural diversity: multivariate analysis. REVS analysis yielded a model with seven predictors significantly related to the extent of stereotyped behaviour (Table 2). Stereotypy showed increasing trends with density of logs on the ground, age and higher among pandas in zoo 3 compared to that of zoo 2, but showed declining trend with number of nests, sociality, tree density and tree height used by pandas (p < 0.05). Together, these seven predictors explained 86% of the variations in the extent of stereotypy among captive pandas in Indian zoos. For behavioural diversity, the REVS analysis yielded a model with eleven predictor variables (also shown in Table 2). Behavioural diversity was lower among individuals in zoo 2 compared to that of zoo 1, during summer than in winter, and decreased with ambient temperature, stereotypy, tree density and tree height used by pandas (p < 0.05). The remaining four predictors in the model were not statistically significant. These eleven factors together explained 79% of the variations in behavioural diversity.

Discussion
Influence of independent variables on the extent of stereotyped behaviour. The overall level of stereotypy we observed was low, suggesting that the pandas in our study were not seriously stressed. The variables that we found to be correlated with stereotypy are consistent with what we know of pandas' natural history. Our study reports that variables like logs on the ground, nest, sociality, zoo, tree density, age and tree height used by pandas are the driving force for stereotypy in captive pandas involved in the study.
Making the captive environment more naturalistic by integrating enrichment into the enclosure seems to be a promising way of alleviating stress and improving both welfare and reintroduction success 41 . It also helps to improve reproductive rate and overall health 39 . Improved health reduces stress and gives greater control over the environment increasing the chances of survival and longevity both in captivity and following release into the wild 5 . It is generally accepted that enrichment of the captive environment increases animals' ability to cope with challenges and positive use of the environment reduces or eliminates aberrant behaviour 23 . Lack of enclosure enrichments and less complex enclosures can cause stereotypy and other atypical behaviours 24 , while providing enrichment increases the frequency of natural behaviours 25 and thereby reduces stress, which in turn decreases stereotypy 27 . But enrichment needs to be appropriate for the species of animal concerned. Abnormal behaviours are often associated with captive conditions that deviate greatly from the species' natural environment. Consistent with this argument we found that though dead and fallen logs on the ground are one of the important characteristics of the panda habitats in the wild [42][43][44][45] , merely providing them in captivity does not ensure the species' welfare: in fact, stereotypy increased with log density in our study subjects. This could be due to the fact that four individuals that showed more stereotypy were housed in the small barren enclosures with no trees but more logs as a part of enrichment. Without those four individuals, the linear relation between stereotypy and log density was not statistically significant. This clearly suggested that merely providing logs in the small enclosures does not maintain welfare.
When animals are housed in enclosures designed to resemble their natural habitat by considering their natural history (provision of vegetation, shelter, pool, etc.), there is a reduction or elimination of abnormal patterns of behaviour such as stereotypies, increased fitness and improved health, all of which may influence reproduction 25,[46][47][48] . For many species, nests, shelter or burrows in enclosures will serve as retreat and hiding places, which are essential to cope with environmental stressors 10 . Gerbils, mice and rabbits have all shown less stereotyped behaviour when retreats are provided 9,49-51 . Such retreats can mitigate the effects of zoo visitors, who can serve as a source of stress for species that rarely interact with humans in the wild. Consistent with these previous results, we found that with provision of nests, the extent of stereotypy decreased in captive pandas. Many species prefer nests both for rearing the young as well as for resting and shelter, and pandas follow this pattern, so providing nests in adequate numbers will supports their natural behaviour as well as provide relief from environmental stressors. Zidar recommends providing one more nest than there are individuals in an enclosure 52 .
Although pandas are an asocial species, our study showed that pandas show more stereotyped behaviour when housed alone than when with another individual or in group. Being a solitary species in the wild might Table 2. REVS model to explore the effect of biological and environmental factors on the intensity of stereotype and level of behavioural diversity among captive pandas in India zoos. Significant values are in bold. *Zoo 1 kept as reference category. **Zoo 3 kept as reference category. www.nature.com/scientificreports/ encourage management to house them singly in captivity, but not every activity and habit of species in the wild can be used in captivity. For example, polar bears are also a solitary species, and it was at one time thought best to manage them alone, but it was found that managing them in a social setting reduces stereotypic pacing behaviour 53 , consistent with this study. Importantly, managers of zoo should note that living in group is greatly influenced by the individuals' compatibility and hence this should be kept in mind while pairing. Similarly, we found that the presence of trees, and greater mean tree height use by pandas, reduced stereotypy. Pandas' preferred high elevation habitat is favourable for taller trees 20 , and Shrestha et al. found that canopy cover was an important factor in habitats for pandas in the wild 54 . In European zoos, pandas spend 90% of their time off the ground 37 . Consistent with these previous findings, our study reveals that more and taller trees support natural behaviours in panda. The Central Zoo Authority (CZA) of India enrichment manual recommends taller tree provision in panda enclosures, and again we provide empirical support for its recommendation.
We found that with increasing age stereotypy increased in pandas. The older the individuals the more time spent in captivity with its associated risks of stereotypic behaviour. The same trend has been observed in other species: for example in captive bears stereotypic behaviour increased with age 55 . In another study Asiatic black bear and sun bear showed more stereotypy with age 56 .
Influence of independent variables on behavioural diversity. As noted in the "Introduction" section, in a species like the panda, high daytime behavioural diversity is not necessarily a positive indication of good welfare. However, our comparison of behavioural diversity with stereotypy showed a negative trend (though not significant), suggesting that low behavioural diversity might be associated with poorer welfare.
Nonetheless, we found some results that suggested that lower diversity might in fact be associated with a more natural lifestyle. Because of the amount of time that wild pandas spend foraging 57 and sleeping or inactive, they cannot show much behavioural diversity, and in our sample of captive individuals, they showed the same trend. For example, behavioural diversity was lower when pandas were provided with more trees in the enclosure. This suggests that when appropriate conditions are maintained in captivity, panda prefer to be inactive during the day, as is consistent with their natural history 57 . As pandas are essentially arboreal mammals, naturally they also spend most of the time inactive (e.g. sleeping) on the trees 57 . Indeed, providing larger trees would promote inactive behaviours and hence lower behaviour diversity in captivity, this captures their natural behaviour. This is consistent with our results where increased tree height used by pandas decreased behavioural diversity.
We found behavioural diversity was greater when there are more logs in the enclosure. In the Yele Reserve in Sichuan, China, Wei et al. found 107 of 185 panda dropping sites (57.8%) on shrub branches, 49 (26.5%) on fallen logs, and only 29 (15.7%) on the forest floor 44 . Droppings were found mostly on elevated structures ranging from 1 to 3 m above the forest floor and occasionally on trees over 12 m. Moreover, microhabitats selected by pandas were also characterized by fallen logs and tree stumps 42,45 . Wei and Zhang mention that to access bamboo leaves easily, pandas usually use some elevated objects, such as shrub branches, fallen logs, or tree stumps to lift their body 43 . Hence, providing tree logs in the vicinity supports their natural behaviour. But at the same time management should keep in mind that merely providing logs in the enclosure would not guarantee species welfare, as discussed in previous section with respect to stereotypy.
Temperature is an important element of microclimate for animals, and influences the activity level of captive animals 10 . When temperature rises, many species show distress in captivity 10 . The red panda inhabits lowtemperature areas 20 , so it is unlikely that higher temperatures would support natural behaviours. We found that with increased temperature behavioural diversity decreased in captive pandas. Similarly, we found that pandas showed higher behavioural diversity in the winter season, where temperatures are low as compared to summer season.
Studies that have tried to relate behavioural diversity and stereotypy in captive animals have varied in their interpretation; many have found significantly inverse relationships between the two 19 . In this study our multivariate model suggested that behavioural diversity is negatively influenced by stereotypy in captive pandas, confirming previous research.
Other factors associated with variations in behavioural diversity are less easy to identify with welfare, positive or negative. Behavioural diversity also decreases with age of pandas and increases with distance to cage mate, number of visitors and quantum of bamboo provided, though these effects were not significant in the REVS model.
Taken together, these results suggest that higher behavioural diversity is not a straightforward indicator of better welfare in all captive animals. The overall non-significant relationship between stereotyped behaviour and diversity we observed could well be the result of a mixture of factors operating in opposite directions. To interpret diversity correctly, it would be helpful to know what level of diversity the species shows in the wild, and such data are rarely available-a limitation of our study as of many others. Although there are dissenting voices 58 , arguably what matters most both in terms of welfare and in terms of potential reintroduction to the wild, is that a captive animal's time budget approximates as closely as possible that of a wild animal. It is not diversity as such that is important, but the behaviours that the animal exhibits.
Differences between zoos. Our study showed that both the extent of stereotyped behaviour and behavioural diversity varied significantly among zoos. However, Zoo 2, an important breeding centre, housed only a female and her two cubs; this may lead to many factors being confounded and is thus a limitation to our study. Captive animals rely on the zoo environment, its routine and husbandry practices to limit their stress levels, and any failure to provide suitable resources will certainly disturb them and lead to distress 10 . Controlling such variables appropriately will help reduce stress among captive animals, and we can rely to some extent on our knowledge of the species' natural history to guide us through this challenge. Our study was able to identify some of the factors www.nature.com/scientificreports/ that are associated with better welfare, but even with these factors taken into account, significant differences among the three zoos remained. These are presumably due to subtler variations in the zoos' environment or management regimes. Since the panda is endemic to high elevations, we considered whether differences between the elevations of the zoos might be relevant, but the biggest differences were between Zoos 1 and 3, which are at essentially the same elevation. In Zoo 1 pandas showed lower stereotypy and higher behavioural diversity then the other two zoos. Again, these differences may be due to subtle differences between the management regimes in the three zoos; possibilities include keepers' attitudes and the zoo's experience in managing pandas. It is notable that Zoo 1 has longer and wider experience in the management of red pandas than the other two zoos, which have joined the captive breeding programme more recently and have fewer animals. Other notable differences were that in Zoo 1, pandas are fed twice a day as compared to the other two zoos where feed is given all at one time (both bamboo and supplementary diet); and that in Zoo 1 the enclosures were of a good size for a small mammal like the red panda, and were well maintained with much natural vegetation. The other two zoos had a large enclosure with poor vegetation (trees and grass), or a small enclosure with a barren floor and no trees at all. Location of the enclosure also needs to be considered: in two of the enclosures at Zoo 3 the sun shone directly on the animals with no shade as such, keeping the temperature higher than would be natural for pandas. Any of these factors could be the reason the pandas performed comparatively well in Zoo 1, and it would be necessary to study a wider (and, therefore, cross-national) sample of zoos holding pandas to identify which of them are the most important.

Conclusions
We conclude that maintenance of adequate number of nests, more and higher trees and log density in the enclosure will reduce stereotypies among captive pandas to a considerable extent. Similarly, encouraging pandas to spend more of the daytime in feeding and resting by providing proper enrichment will support their natural behaviour. It is notable that these factors are also important habitat attributes for pandas in the wild 54 . Our study is important as the first to provide detailed empirical study of the welfare of captive pandas, and factors influencing it in the Indian subcontinent, where a natural population exists. A limitation of this study is that there are no comparably detailed studies among panda populations in captivity as well as from the wild around the world. For captive breeding purposes, individuals are exchanged among almost every participating breeding centre in the world (Darjeeling Zoo Management Plan), and in the process individuals with welfare problems will also be transferred. Ultimately this can reduce the success of the expensive global breeding programme. The Indian captive breeding programme holds especial potential because its sites are within the species' natural range. Its primary responsibility is to maintain the species in the wild through reintroduction. At several points our study provides empirical support for the recommendations of the Indian CZA, and as captive red panda populations experience similar welfare issues among breeding centres around the world, the findings from the present study could help to refine the captive protocol to enhance the welfare of pandas at captive breeding centres in India and elsewhere in the world. In brief, we recommend.
• Enough nests (at least one more than the number of animals in the exhibit) should be provided in the enclosure, as a retreat, especially during the breeding season, where rearing young ones is an essential part of the species' survival. • To support species-natural behaviour, and replicate wild conditions, panda enclosures should be enriched with an adequate number of tall trees and multiple tree logs or stumps in a novel, yet complex way with a number of different pathways for simultaneous use by housed individuals. • Feed should be provided at different time intervals at least twice a day.
Although we did not find a significant relation between behavioural diversity and stereotypy, our multivariate analysis has suggested that stereotypy would affect behavioural diversity inversely in captive pandas. We suggest studying in detail the relationship between behavioural diversity and stereotypy with respect to captive pandas.

Method
Study sites. The study was carried out at three different Indian zoos that are part of captive breeding programme (Table 3). More details about the zoo, subjects and their husbandry practices are given in Supplementary material 1.   Table 4). We also took into account the zoo literature on red pandas, and discussions with zookeepers at the study sites, while preparing the ethogram. To confirm the credibility of the ethogram we compared it with international literature on panda behaviour 37 . We divided the behaviour of the panda into active, inactive and stereotyped behaviours, using the ethogram (Table 4). During observation, when the subject went out of the observer's sight, we recorded it as ' out of sight' (mean 0.64%, range 0-5%). We calculated the percentage of time spent on stereotyped behaviour for individual animals to compare with the independent factors.
Behavioural diversity calculation. We measured behavioural diversity for each individual using the Shannon-Wiener diversity index H 18 using the formula:   Table 5. In the case of behavioural diversity, we included the extent of stereotypy as an additional independent variable. We have included season as one of the factors, owing to significant difference in ambient temperature between winter (as low as zero) and summer (as high as 25 °C), with which the pandas alter their behaviours. Accounting season as one of the variables reflects this behaviour adaptation if any. The season considered in the study represent the actual observation time and its effect on the behaviour.
Statistical analysis. Statistical analysis was carried out using SPSS for Windows (Version 16) and R (Version 3.4.1). Before testing for differences in the levels of stereotypy and behavioural diversity due to various biological and environmental variables, we checked the data for normality. The data for the dependent variables (stereotypy and behavioural diversity) were not normally distributed, so we performed non-parametric tests. For the behavioural diversity data set we excluded two outlier data points with respect to the tree density variable. The univariate dependence of stereotypy and behavioural diversity on biological and environmental factors was tested using the Mann-Whitney U test (for factors with two categories) and Kruskal Wallis H test (for factors with three categories). For the season category (winter and summer) we coded subjects as 0 and 1; 0 for winter observation and 1 for summer season. The first-order effects of the continuous variables were checked using simple linear regression; in particular, to check the effect of stereotypy on behavioural diversity, we treated the former as independent variable and latter as dependent variable in linear regression. To quantify and assign empirical support to the simultaneous effects of each of the biological and environmental factors, we employed H = − i p i ln p i